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SUMMARY 

A mixed analysis-inverse procedure based on the full potential equation 
in conservation form has been developed to recontour a given base wing to pro- 
duce a prescribed favorable pressure distribution. The method incorporates a 
novel density linearization scheme in applying the pressure boundary condition 
in terms of the velocity potential. The FL030 finite volume analysis code has 
been modified to include the inverse option. The new surface shape information, 
associated with the modified pressure boundary condition, is calculated at a 
constant span station based on a mass flux integration. The inverse method is 
shown to recover the original shape when the analysis pressure is not altered. 
Inverse calculations for weakening of a strong shock system and for a laminar 
flow control (LFC) pressure distribution are presented. Two methods for trail- 
ing edge closure model are proposed for further study. 

INTRODUCTION 

Currently, the aircraft industry is in need of quick turnaround methods 
to develop energy efficient transonic configurations with optimal aerodynamic 
characteristics. Development of computational transonic methods over the last 
decade has significantly contributed towards fulfilling this need by aiding 
the design of efficient transonic airfoil sections and wing surfaces. Although 
computational models have been primarily developed to treat the direct problem 
of determining the load characteristics of a prescribed shape, the inverse 
problem associated with determining the required recontouring of a given wing 
to provide a preassigned favorable loading is becoming increasingly important 
to eliminate much of the cut-and-try approach to geometry definition. 

Inverse methods based on the transonic small disturbance theory^^^ in 
two^^^ and three^^"^^ dimensions and full potential models in two^^"®^ and 







three^^^ dimensions have been developed with some restrictions or other. The 
small disturbance method^ ' provides geometric versatility in designing fairly 
arbitrary geometries. However, the limitation of the method involves the 
breakdown of the theory for large flow deflections, especially near the 
leading edge. The existing full potential inverse method^^^ that can handle 
the design of shocked flows is based on the nonconservative form of the full 
potential equation and uses the FL022 analysis code^®^. It is essential that 
the finite-difference approximation to the full potential equation be cast in 
conservation form to satisfy certain jump conditions' ' across the shock 
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system. The nonconservative procedures' introduce mass sources at shock 
waves, and the strength of these sources depends on the local grid spacing, a 
non-physical consideration. Erroneous shock solutions could thus result in 
improper geometry definition while using inverse methods based on nonconserv- 
ative formulation. 

Other inverse methods such as the ones based on the "fictitious gas" 

approach^^®"^^^ are oriented toward achieving shockless designs. Such a 

restriction may be too severe from the standpoint of aerodynamic efficiency, 

since some wave drag may be necessary for the production of a good lift-to- 

drag ratio. Of equal significance is the fact that a shockless wing could 

experience radical trim changes associated with sudden generation of large 

aerodynamic center shifts produced by shocks at slightly off-design conditions. 

In general, inverse methods provide a valuable alternative to optimization 
(131 

methods' ‘ which can provide shapes that optimize certain aerodynamic quan- 
tities but require excessive computer time for any realistic wing modification. 

The present report deals with the development of an inverse method based 
on the fully conservative form of the full potential equation to address some 
of the limitations of the existing methods. The currently available FL030 
finite volume full potential analysis code for wing-body combinations is 
modified to include the inverse option. The crux of the inverse problem is 
the incorporation of the prescribed pressure as a boundary condition on a 
surface yet to be determined as part of the solution procedure. A density 
linearization scheme is introduced in this report in applying the pressure 
boundary condition in terms of the velocity potential. Initially, the pressure 
boundary condition in terms of a Dirichlet problem is applied at the original 
shape location. After every n iterations (n 5), the new shape information is 
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obtained at every span station using a mass flux integration procedure. 
Application of this inverse procedure to v/eaken the shock system of a typical 
transonic wing is illustrated. Another example of a wing design for laminar 
flow control pressure distribution is also demonstrated. The inverse method is 
reasonably inexpensive (35 to 45 minutes of CDC 7600 time for an analysis- 
inverse calculation) to use for wing modification requirements. The inverse 
program is also operational at the NASA-Langley Research Center using the 
CYBER 203 computing system. 

At present, the currently developed inverse code is only a research tool 
and requires much more work to understand the constraints to be imposed on the 
specified pressure to achieve physically realistic looking shapes with closed 
trailing edges and also the relationship between the freestream Mach number and 
the specified pressure. 


FORMULATION 

The conservative form of the full potential equation in a general coordi- 
nate system ?,n»? can be written as shown in Eq. (1) below. (This report uses 
(x,y,z) -*■ (?,Ti,?) as notation for the transformation. The use of (x,y,z) -+■ 
(?>n>c) is also common in the literature.) 

where U, V, and W are the contravariant velocity components, p is the density, 
and J is the Jacobian of the transformation that relates the general coordi- 
nates ?,n,5 to the Cartesian system x,y,z. Introducing the following notation 
for convenience 

= U , Ug = V , U 3 = W 

x^ = X , Xg = y , X 3 = z 

= C , Xg = n . X 3 = E 

the contravariant velocities are given in terms of the velocity potential <|) by 
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U. = Z i=l,2,3 

' j=l ^j 




3X. 3X. 
_L _J. 
3X^ 3X^ 


i=1.2,3 

j=l,2,3 


The Jacobian of the transformation J is represented by 



The density p is computed from the isentropic formula 


p = [. 


^ M“{q= 


')] 




( 2 ) 


(3) 


(4) 


where the total velocity q is obtained from the relation 
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An analysis problem is one in which the Eq. (1) is solved to produce the 
flow field over a given geometry by imposing the usual surface tangency 
boundary condition = 0 (n is normal to the body surface) on the exact surface 
location. If n is the coordinate leading out of the surface, then the surface 
tangency condition reduces to the simple form in terms of the contravariant 
velocity V 

V = 0 (6) 

on the surface. After Eqs. (1) and (6) are solved together, the resulting 
pressure distribution over the surface is computed from 
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(7) 


C 


P 


(p^-1) . 


An inverse problem is one in which the Eq. (1) is solved subject to a 
prescribed pressure distribution (Cp specified) and the resulting body shape 
that satisfies the surface tangency condition Eq. (6) is sought. 

Usually, for easy handling of the boundary condition, a body fitted 
coordinate system is chosen for Unlike the analysis boundary condition 

(V = 0), the incorporation of the inverse boundary condition in terms of a 
prescribed Cp (Eq. (7)) is considerably more difficult because the velocity 
potential <}) appears nonli nearly through the term in Eq. (7). In order to 
aid in the application of the inverse boundary condition, first the density p 
appearing in the Cp relation is linearized as follows. 

Density Linearization 
From Eqs. (4) and (7), we can write 



It can be seen from the above nonlinear relationship that from a given Cp 
distribution extracting the information in terms of the velocity potential <j) 
would involve some type of a linearization. Denoting the current iteration 
cycle by (n+1) and the previous one by n, the variation in density due to 
variation in <p can be expressed as 

n+1 - n . . 

p = p + Ap (9) 

where p*^ = p(cj)'^), Ap = p(c|)'^ + Acj)) - p(<|)'’) and Acf) = - cj)*^). Rewriting 

Eq. (8) fully in terms of (j), p((ji) can be expressed as 
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( 10 ) 



Substituting ((j) + A(j)) in the place of 4» in Eq. (10) and using binomial expan- 
sion, an expression for p((j) + Act)) can be written as 


^ p" - (p")''V[u" A . V" + W" ^1 A* . (11) 

The derivation of Eq. (11) is given in Appendix A. 

xU 

While operating at the (n+1)^ iteration cycle, all the quantities appear- 
ing at the n^^ level are known and Eq. (11) can now be used to get an estimate 
for A(j) at the body surface for a given pressure distribution. Since we require 
V = 0 at the body, the given Cp can be expressed as 



differential specified 


operator 

In the inverse problem Eq. (12) will be discretized to get an estimate 
for A({) = (<f)'^^^ - (j)*^) which in turn will be used as a Dirichlet boundary 
condition while solving Eq. (1). 


Implementation of Boundary Conditions 

When Eq. (1) is discretized and written in terms of Acfi using Jameson's^^^^ 
pseudo-time concept, at any point (i,j,k) it will appear in tridiagonal form as 




(13) 


where TM, T, and TP are the coefficients of the tridiagonal system with built 
in artificial viscosity for handling mixed elliptic-hyperbolic flows and R is 
the finite-difference operator to be satisfied and is evaluated using 
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values of c(> from the previous iteration and values of (J> which have already 
been updated on the current iteration. Referring to Fig. 1, at any 
boundary point (• symbol) the evaluation of TM, T, TP, and R would require 
velocity potential information at the dummy points (□ symbol) that are intro- 
duced inside the body surface. Boundary conditions on the surface play a role 
in eliminating this dummy point information. 

Analysis Problem 

The analysis problem imposes V = 0 at all body points by simply reflecting 
all the various flux quantities across the surface. Referring to Fig. 1, this 
is done by setting 


O ~ FIELD POINTS 
• ~ BOUNDARY POINT 
□ ~ DUMMY POINT 


SECONDARY 




Fig. 1. Boundary cell distribution 
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(" 'A,u ' ' (' 51...I 


(•si.,, 


(•51.,, ■(•SI,. 


Equation (14) would automatically set V*^ = 0 while forming TM, T, and TP, 
but doesn't rigorously satisfy =0 which is the actual boundary condition 
to be imposed. This can probably be achieved if corresponding to 

the dummy point can be replaced in terms of information on the surface and 
above the surface appropriately. In the present method (A4 >).j is simply 

set to zero while solving for the body point. 


Inverse Problem 

Referring to Fig. 2, when Eq. (13) is written at one point above the body 
surface (point A at i,j-l,k), it involves (A4 >)j _• ^ appearing at the body 
point. In the inverse problem, the value for (A(ji),. . f, at the body point is 
first computed from the prescribed pressure distribution using the density 
linearization procedure given by Eq. (12), in the following way. Referring 
to Fig. 3, the pressure coefficient is prescribed at the center {* symbol) of 

SC82-18227 


Fig. 2. 

Grid point notation 
for the inverse 
procedure 
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each primary cell face coinciding with the body surface. First consider the 
lower surface where along the direction of sweep the index i increases. The 
discretized form of Eq. (12) can be written as (at point P in Fig. 3) 



Since the direction of sweep is along increasing k-index in the span direction 
and increasing i-index in the streamwise direction at the lower surface, the 
quantities (Acj))^, (A4 ))q, and (A<t))j are known and the unknown to be computed 
from Eq. (15) is (A<{))<; = (Ac()). . . . This is required while solving Eq. (13) 
at point A in Fig. 2. On the upper surface where the i-index is decreasing 
along the direction of sweep, Cp prescribed at (i+^^.k-^s) is used to compute 
(M).- j 1 ,. For example (in Fig. 3) the pressure coefficient at point N and 
(A(j))^, (A(J))q, and (A(|))j^ will be used to compute (Atj))^^ = (A({))^. j in a manner 
similar to the Eq. (15) for the lower surface. While solving Eq. (13) at 
point A in Fig. 2, the quantity TP(A4>). ^ i, is known from the above procedure 
and is lumped into the right-hand side residual term and Eq. (1) is solved 
only up to one point above the body surface. Thus, the inverse problem uses 
a Dirichlet boundary condition. 

New Shape Information 

Initially, the pressure boundary condition is applied at the original 
shape location. After every n inverse relaxation cycles (n 'b 5 to 10), the 
new shape information is obtained by using a mass flux integration procedure 
as follows. 
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Referring to Fig. 2, point B is on the old surface where the specified 
pressure condition, in terms of was imposed as a Dirichlet boundary 

condition. After a vertical line relaxation is completed, the finite differenced 
form of Eq. (1) given by Eq. (13) is solved at point B, using (A<{))g and (A(|))y^ 
now available. The dummy point value of A4) in Eq. (13)) is set to zero, 

just as in the analysis problem. The right hand side R in Eq. (13) at point B 
can be represented as R = R ||p , (p analysis calculation 

jj is set equal to ”|p j) • 3*^ inverse problem, where the new 

shape information is sought, the flux value |p will not be equal to -|p . 

By accepting the value for ^p as it exists at point C, solution to Eq. (13) 
at point B will yield a value for the flux |p . The modified flux information 
at the old surface point B is taken to be ^p = -^j^p + ^p |. Again, 


this will not be zero for an inverse calculation. Once the modified flux 
information is known at the old surface points, the new shape information 
can be obtained. Let the dashed line in Fig. 4 represent the modified 
new shape. The surface transpiration at i-1 grid point is denoted by (dri)^-_i> 
and at point B by (dn)^-. Balancing the mass flux between the old shape 
(solid line) and the new shape (dashed line), the following relationship is 
obtained (neglecting the effect of the spanwise variation) 



(• 5 ),., 


( dn )^_1 I 



- Ci_l)= 0 . (16) 


Equation (16) assumes that V is zero along the dashed line (boundary condition 
for a solid surface). The only unknown in Eq. (16) is (dn)^. Usually, the 
nose shape is prescribed, and the starting value of (dri)^._^ is zero at the 
point of transition from analysis to inverse. Once (dn)^- is known, the new 
values of x and y at point B" are computed as follov/s: 




where (x^)^ and (yj^ 


are obtained by three-point one-sided differentiation. 


RESULTS 

The finite volume FL030 code^^^’^®^ is an analysis code based on the full 
potential equation in conservation form and has the capability to handle wing- 
body combinations. The inverse procedure presented in this paper is also based 
on the full potential equation in conservation form and the FL030 analysis 
program is found to be a good choice to incorporate the inverse logic. One 
advantage of using the FL030 program is that it requires only a local descrip- 
tion of the coordinate mapping to a body-fitted system and essentially 
decouples the solution process from the generation of the grid network. As 
a result, during the inverse calculation as shape changes take place, this 
method requires grid adjustments only to local cells adjacent to the wing 
rather than having to change the entire grid distribution at the end of each 
relaxation cycle. 

To test the inverse concept, first an analysis calculation v«s performed 
using a typical transonic wing geometry definition as shown in Fig. 5, at 
= 0.86 and freestream angle of attack of 4.68°. After a sequence of crude- 
medium-fine grid calculations (approximately 30 minutes of computer time on 
the CDC 7600 machine using a 161 x 27 x 35 fine grid), the analysis calculation 
was reasonably converged. The resulting pressure distributions on the upper 
surface at discrete span stations are shov/n in Fig. 5. The presence of a shock 
system is evident and the strength of the shock gradually increases from the 

wing root reaching a peak strength around 85^ span. As a verification for the 

correctness of the inverse procedure, the analysis pressure of Fig. 5 is kept 
unaltered and specified as input pressure for the inverse calculation. After 
20 inverse cycles, the resulting shape information is provided in a tabular 
form from the computer output in Table 1. It has eight columns. Explanations 
for Columns to are given below. 

Column Value of x/c at that span station. 

Column Value of x of the surface grid point. 

Column Q): Value of Y of the grid point on the original shape. 
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Fig. 5. Typical transonic wing showing presence of a shock system at 
= 0.86, a = 4.86° 

00 

Column (^: Value of z of the surface grid point. 

Column Cp at node point on the surface (i,j,k). 

Column d): Cp at half node point (prescribed). 

Column d); Value of Y of the grid point on the new shape. 

Column d): Index i in the ? direction. 

In Table 1, the Cp in column (J) is the same as the analysis calculation 
of Fig. 5. The resulting shape information in column (?) closely 
duplicates the original shape given in column (J). 
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Table 1. Recovery of original shape for unaltered analysis pressure _ 
specificatlont “ 0.86, o = 4,86“. (Explanations of columns 
through (£) are on pages 12 and 13.) 

SECTION CHARACTERISTICS 


PCT SEKISPAN CL CD CM 

_.5614^ -3348 ^ .^36 ^ -.1695 

© ® ® ® AJ) ® ®“ 

1.000 559.301 118.567 127.750 .rTfc .185 llciT620 29 

.967 555.793 116.776 127.711 .141 .167 118.931 30 

. 9 34 5 52.35 - fc ^ l -9r i -3 ^ - ie - 7» - 6 »6 .109 119 . 1 9 1 31- 

.902 548.916 119.491 127.660 -.073 -.030 119.543 32 

.870 545.538 119.736 127.647 -.128 -.132 119. 7BB 33 

.83 9 542.1 -9 6 -1 1 -9^5 1 - 7^^7. - 6 4 1 - r - 156 - .1 4 1 - 11 - 9. - 969 34- 

.808 538.894 120.095 127.635 -.145 -.145 120.149 35 

,777 535.632 120,251 127.631 -.158 -.159 12C.304 36 

.7 4 7 532. 4 10 - 12 -e y 4 {Ki H t2 - 7r62 - 0 r^-89 rl 6r 9 - 12 - & .-45g 3^ 

.717 529.229 120.553 127.601 -.290 -.209 120.607 38 

.687 526.094 120.661 127, 5B7 -.458 -.389 120.715 39 

. 6 5 8 -5 2 -3 . 0 0 4 120 . - 7 5- 9 ~ 1 2 ~ 7 . -5 75 - . -5 65 r5z r 2 - i 2- l -rB^tt 46- 

,630 519,961 120.836 127,565 -.617 -.594 120. B90 41 

.601 516.967 120.881 127.560 -.630 -.618 120.936 42 

— . -5- 7 -4 -5 1 t4- » 02 4 120.898 127.55 9 - . - 622 .619 - 12 - 0^ 95 2 43 

.546 511.134 120.876 127.563 -.608 -.605 120.930 44 

.520 508.297 120.835 127.570 -.599 -.591 120.999 45 

— . 49 3 5 0 3 . 3 1 H^ -1 2 3 *776 -12 - 7. 5 79 - . S OL ■ y 59- fi - 130r E> 2 9 46- 

.468 502.768 120.687 127.590 -.570 -.569 120.739 47 

.443 500.119 120.571 127.603 -.553 -.549 120.622 48 

■ 4- 1 8- -4 9 3 . 5 0 9 12 0.4 45 - i 2 - 7 .6 1 - 7 - r54 - ? r 5 3-9 - 1 2 - 0W t-95 49- 

,394 494.960 120.301 127.634 -.536 -.530 120.350 50 

.371 492.472 120.146 127.652 -.535 -.525 120.194 51 

— .348 49 0.049 119.979 - 12 - 7.671 - .538 r 5 23 — l -e- O rfrgfi 5P- 

.326 487.690 119.802 127.692 -.543 -.532 119.349 53 

.304 485.398 119.612 127.714 -.544 -.538 119.656 54 

— . 2« 3 - 4 - 83 - . 1 - 7 4 - 1 - 1 '9 . ^ 0^ - i2^» - ^3 -6 r5*9 . 5 5 5 - il - 9r 45 0 55- 

.263 481.020 119.186 127.764 -.537 -.528 119.228 56 

.243 478.936 118.965 127.799 -.539 -.532 119.005 57 

— .22 4 4T6 - . - 927 118.723 127.91^ - . 5 3 5 .554 - ll - r» - 7&2 55- 

.206 474.993 116.479 127.847 -.530 -.525 118.516 59 

.189 473.135 118.227 127.676 -.525 -.522 118.263 60 

— .172 471 - *35 4 11 7. - 9 5 5 ^g - 7r 9 (H? r5rt5 » - 5 15-1 1 - 7 .- 9 9 C 64- 

.156 469.651 117.694 127.937 -.503 -.503 117.726 62 

.140 468,027 117.418 127.968 -.491 -.490 117.449 63 

— » 126 - 4 66 - . 4 e - 3 -- l - l - 7. - l r 45 - ig7r 9 9 € - < r48- 7 r4 8- 2 — 1 -4- 7 -r 1 -7 1 64- 

.112 465.019 116.869 128.026 -.493 -.493 116.695 65 

.099 463.635 116.599 128.054 -.504 -.491 116.622 66 

— > 6 8 7 - - 4- 62 - * -5- 3 5- - M - 6 . -5 2 -4 - ^28 »0D - 2 r568 r5 6 5 -i -4- 0 i- 5 -4 6 67- 

.075 461.119 116.044 128.109 -.507 -.502 116.064 68 

,065 459.985 115.767 126.136 -.507 -.504 115.795 69 
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Table 1 (Continued) 

.046 457.97? 115.212 126.187 -.515 -.510 115.226 71 

.038 457.093 114.931 128.212 -.523 -.514 114.942 72 

■ *i 56 - »29 5- li -«( » 6rb C - 1281235 - . 5 3 9 - »5£ 9 1 - 1 4.66 7 ?d- 

.023 455.581 114.377 128.259 -.533 -.546 114.394 74 

.017 454.953 114.084 128.263 -.510 -.517 114.088 75 

- -.- 01g -4 54 - .3» ^- ll -- 3.7 ^- 3 12 6 . - 306 -.48 2 .503 1 - 13 . 7 9 4 ?6r 

.008 453.933 113.476 128.330 -.395 -.459 113.476 77 

.004 453.548 113.135 128.355 -.235 -.331 113.139 76 

- » - a 0 2 - 4 5 5.2 7 - 7 - 1 1 - 2 « 76 r 1 2 0 « -3 9 2 .131 1 - 12.7 6 1 IS- 

.000 453.142 112.352 128.410 .264 .147 112.352 8C 

0. 453.109 111.944 126.437 .457 .375 111.544 61 

" t OO t —4 53 - . - l 7 - i-l M - « 54 2 12 6 . 4 62 .-56 t8 . 526 111 . 542 8£- 

.002 453.325 111.146 128.485 .593 .594 111.14B 83 

.004 453.576 110.766 128.507 .553 .566 110.766 84 

*008 453* 9 22 11 0.3 9 7 128*526 r454 *51 5 110*397 frir 

.012 454.373 110.055 128.544 .340 .401 110.061 86 

.017 454.913 109.733 128.559 .230 .292 109.750 87 

. 023 45 5-.54-5 -1 0 9. 4 - 31 12 8. 573 .-t2-7 . 1 6 5 109 .4 5 6 86^ 

.030 456.276 109.170 128.585 .045 .069 109.196 89 

.038 457.095 106.534 126.595 -.012 .021 106.964 90 

• .046 -- 4 5 e.0DD 100.723 120. 6 0 4 -.0 5 0 .02 7 10 6 . 7 5 6 91- 

.055 458.990 108.531 126.611 -.C79 -.057 106.566 92 

.065 460.062 108.349 128.616 -.108 -.085 108.387 95 

■ ■.- 07 6- 461 . 219 10 6 .1 8 7 126 .6 24 —- 13 0 -. 115 106 . 22 6 94- 

.088 462.456 108.041 126.630 -.144 -.131 106. 0B4 95 

.100 463.779 107.907 126.634 -.159 -.143 107.953 96 

--. 11 4 46 5 *102 107. 78C 126.639 - »17 9 - .1 6 1 107.8 2B 97- 

.128 466.662 107.668 128.642 -.197 -.162 107.718 98 

.142 468.224 107.574 126.646 -.210 -.197 107.627 99 

-. 15 6- 4 69.66 4 197.4 9 2 • 1 - 2 6.648 223 — 2 0-9- 1 - 07 . 5 4 7 3r0-0- 

.174 471.582 107.426 128.650 -.235 -.222 107.433 101 

.191 473.376 107.373 126.552 -.247 -.233 107.432 102 

- .200 4 7 5 .2 45 107.3 3 6 - 120 . 653 —- 25 6 - .? 4- 6 1-0 7 . 3 9 6 1-8^ 

.227 477.186 107.317 128.554 -.268 -.252 107.379 104 

.246 479.201 107.299 126.656 -.289 -.270 107.362 105 

-r 26 5-48 1 - .2 86 10 - 7. 5 - 10 13 6 .6 57 - .505 . 2 9-2 -1 - 0^ . 573 

.286 483.441 1C7.338 128.658 -.317 -.305 107,402 107 

.307 485.665 107.381 128.658 -.328 -.319 107.446 108 

, 52 8 -48 7 . 95 7 ■ 1 07 . 4 4 6 - 1 - 2 8.6 57 -. 541 — - 5 29 -1 - 0 - 7 . 5 - 11 17-9- 

.350 490.314 107.528 128.657 -.357 -.346 107.594 110 

.373 492.736 107.634 126.655 -.371 -.364 107.700 111 

■ .3 9 7 495 . 221 107 . 7 6 5 12 8.6 53 -.-524 -.- 3 52 ■ l - C - 7 .g- 31 14r2- 

.421 497.767 107.922 128.649 -.394 -.394 107.988 113 

.445 500.373 106.115 126.644 -.399 -.403 108.182 114 

. 4 70 5D 3- .03 ^- l & C* 3 4 2 - 12PrC37 - .39 8 -. 4G - 7 1 - 0F.409 

.496 505.761 106.611 126.627 -.363 -.407 108.679 116 

.522 508.540 108.932 126.613 -.350 -.376 105.002 117 

- .54 9 51 1-.- 37 - 2- 169 . 30 4 126 , 5 9 6 .■ 5 5- 7 1 - 09 . 374 irirh- 

.576 514.257 109.737 128.576 -.254 -.239 109. B09 119 

.603 517.192 lie. 238 128.550 -.185 -.225 110.312 120 
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Table 1 (Concluded) 


• 66 0 
.689 

523.205 

526.281 

111.410 

112.060 

128.484 

128.442 

4 ^4 ft T 

-.031 

.045 

4 4 1 

-.073 

.006 
ft 7 Q 

111.488 

112.140 

4 4ft ft ft 1 

122 

123 

4 ft /i 

4 r X c 

.748 

.778 

DC 7 4 ■! 1) U 

532 .563 
535.767 

c. •» ri ft 4 ▼ 

X X £ 4 t c X 

113.384 

114.030 

4 4 i. ^ C ft 

X L • d ./ f 

128.348 

128.289 

< *1 Cl O »» 4 

• X XO 

.173 

.227 

'"■'“•"U’f 

.142 

.199 

»»C ») 

113. 46S 
114.116 

^ 1 /i ^ t c 

X (! *1 

125 

126 

< »i-» . . 

VU U V 

• 643 
.871 

OO J ■ U X o 

542.298 

545.624 

115.224 

115.719 

4 4 r 4 » n 

128.154 

128.091 

^ r\ f. ft » 

.323 

.362 

.301 

.345 

1 ft 

115.314 

115.311 

4 4/- ft ft 1 

128 

129 

' 9*7v O 

.935 
.967 
1-4-09 3- 

o •707 

552.393 

555.836 

559.301 

X Xu • X ^ V 

116. 396 
116.477 
116*-317 

X c. (^4 V U u u 

127.989 

127.963 

.401 

.365 

• W / y 

.403 

.400 



X X tz 

116.490 
116.570 
■116.4 0e ■ 

1 

131 

132 

-1-33 


Next, the analysis pressure distribution with a strong shock system was 
modified on the upper surface from 50% to 95% span in such a way that the shock 
strength is considerably reduced. The modified pressures were then used as an 
input to the inverse code and the inverse calculations were started from the 
converged analysis results. After 50 fine grid cycle inverse calculations 
(15 minutes of CPU time), the residual and the maximum change in the velocity 
potential were of the same order as the converged analysis calculations. A 
sample output of this inverse calculation at 70% span station is shown in 
tabular form in Table 2. The shape differences between the original shape and 
the modified shape can be seen by comparing column (J) and column ©. 

Figure 6 shows the same results in graphical form. The openness of the trail- 
ing edge for the modified shape is much smaller than the original shape. 

One other design problem reported here is a laminar flow control wing 
design. The objective here is to start with the base wing geometry shown in 
Figs. 7 and 8, and then modify the airfoil sections in the test strip shown 
in Fig. 7 to produce a laminar flow control pressure distribution shown in 
Fig. 9 at the midspan region. This is a very difficult design problem because 
the prescribed pressure is considerably different from the one produced by the 
base wing geometry. Like the previous example, the inverse calculation with 
the specified LFC pressure at 50% span was started from the analysis calcula- 
tion. After 50 design cycles (the modified shape was computed at the end of 
every 5 cycles and updated), the resulting modified shapes at two different 
span stations to provide the LFC pressure of Fig. 9 at midspan, are shown in 
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Table 2. Computer output at 70% span indicating the old and nev/ shape for an 

inverse calculation to weaken the shock strength, = 0.86, a = 4.86“ 


SECTION CHARACTERISTICS 


PCT semispan 
-.69 32 , 

CL 

.3414 

CD 

-.0034 


CM 

-.1724 




© 

© 

& 

CD 

(b 

(D 


566.782 

116T347 

157.993 

• T9D 


114^951 

29 

.967 

56 3 . 8 6 fc 

116.449 

157.963 

.163 

.132 

114. 9B5 

30 

.93^ 

560.993 

116.653 

157.923 


• 14 5 

115.169 

31 

.902 

558.148 

11 6.868 

157. BSl 


• 04 D 

115.327 

32 

.870 

555.338 

1 1 7.030 

157.865 

KHbeiH 



35 

.839 

552.563 

117.159 

157.665 

-.149 


115.381 

34 

.8 08 

549 .820 

117.292 

157. B64 

-.203 

-.176 

115.465 

35 

.777 

54 7.11 1 

117.405 

157. P66 

-.263 

-.230 

115.591 

36 

.747 

544.436 

117.515 

157. B67 

-.323 

-.274 

115.814 

37 

.717 

541.797 

117.620 

157.869 

-.390 

-.337 

116.139 

38 

.687 

bl^lAldsI 

imfsTi 

157.673 


-.424 

116.533 

39 

.658 

536.632 

11 7.752 

157. B79 

-.444 

-.447 

1 1 6.996 

40 

.629 

534.107 

117.799 

157.396 

-.493 

-.464 

117.396 

41 

.601 

531.622 


157. B96 

-.572 

-.54 2 

11 7.613 

42 

.573 

529.179 

117.790 

157.889 

-.619 

-.609 

117.700 

43 

.546 

526.779 

117.742 

157. B94 

-.626 

-.625 

117.704 

44 

.520 

524.424 

117.677 

157.901 

-.627 

-.620 

117.664 

45 

.493 

522.115 

117.596 

157.909 

-.624 

-.624 

117. 5B6 

46 

.466 

519.854 

11 7.48 3 

157.921 

-.609 

-.612 

117.466 

47 

.44 3 

517.640 

117. 359 

157.933 

-.595 

-.593 

117.331 

48 

.416 

515.477 

117.221 

157.946 

-.587 

-.53 6 


49 

.394 

513.364 

117.068 

157.961 

-.581 

-.577 

11 7.014 

50 

.371 

511.302 

116.907 

157.976 

-.578 

-.573 

116. B40 

51 

.348 

509.294 

116.737 

157.992 

-.579 

-.572 

116.657 

52 

.326 


116.559 


-.584 

-.576 

116.464 

53 


505.44 1 



-.582 

-.530 

116.265 

54 

.28 3 

503.599 

116.178 

158. C42 

-.572 

-.573 

116.065 

55 

.263 


115.967 

158.056 

-.568 

-.564 

115. B46 

56 

.243 


115.758 

156.074 

-.572 

- . 5 6i. 

115.628 

57 

.225 

496.419 

115.540 

156.090 

— . 566 

-.569 

115.404 

58 

.206 

496.814 

115.312 

158.107 

-.556 

-.556 

115.174 

59 

.189 

495.270 

1 15. 084_ 

156. 124 

-.551 

-.551 

114.947 

60 

.172 

493.791 

114.848 

158.141 

-.541 

-.544 

114.713 

61 

• 156 

492.376 

114.610 

158.158 

-.529 

-.531 

114. 4B0 

62 

.141 

491.026 114.367 


-.520 

-.51 B 

114.245 

63. 

.126 

489.742 

114.127 

156.192 

-.517 

-.513 

114. OlE 

64 

.112 

488.523 

113.689 

158.206 

-.520 

-.510 

113. 79B 

65 

.099 

487.372 

113.655 

158.223 

-.531 

-.513 

1 1 3.533 

66 

.087 

486.289 

113.417 

156.239 

-.535 

-.532 

113.366 

67 

.075 

485.276 

113.175 

158.254 

-.528 

-.527 

113.146 

68 

.065 

484.331 

112.935 

156.269 

-.520 

-.523 

112.931 
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Table 2 (Continued) 


.055 483.458 11?.'651 158.284 -.513 -.513 112.710 70 

.046 482.652 112.452 158.299 -.507 -.511 112.489 71 


.030 

481.25 0 

111.964 

158.327 

-.503 

-.505 

112.014 

73 

.023 

480.654 

111.714 

158.341 

-.516 

-.511 

111.749 

74 

'. 0 1 7 

480 . 129 . 

111.453 

158.356 

-.503 

-.527 

111.468 

75 

.012 

479.669 

111.188 

158.370 

-.442 

-.433 

111.192 

76 

,007 

479.282 


158.385 

-.324 

-.415 

110.904 

77 

. 0 S 14 

478.968 

_ U 0 . 601 _ 


-. 14 * 

-.260 lie . 601 

78 . 


478.753 


158.417 


-.053 

110.266 

79 


478.648 

109.914 

158.434 

.324 

.204 

109.914 

80 


478.625 

1 09.568 

158.449 

.489 

.309 

109,568 

£1 

.001 

478.68 1 

109.229 

156.464 

.577 

.539 

109.229 

82 

.002 

478.814 

108.899 

158.478 

.583 

.591 

108 . B 99 

63 

.005 

479.030 

_ 10 e . 5 SJ _ 

158.491 

. .515 

,562 108,583 

84 

.008 

479.326 

108.284 

158.502 

.395 

.469 

108,284 

85 

.012 

479.712 

108. 018 

158.512 

.266 

. 34 2 

108.019 

86 

.017 

480.165 

107.763 

158.522 

.143 

.213 

1 07.774 

87 

.023 

480.694 

107.529 

158.530 

.049 

.110 

107.553 

68 

.030 

481.303 

107.329 

158.537 

-.023 

.022 

107.371 

69 

.038 

481.984 

107.151 


-.072 

-. 03 ° 

107.203 

90 

.047 

482.736 

106 . 992 

156.548 

-.104 

-. 079 

107.055 

91 

.056 

48 3.55 7 

106.848 

158.553 

-.131 

-.106 

106.921 

92 

.066 

484.446 

106.714 

158.557 

-.157 

-.154 

106.795 

93 

.077 

485.405 

106.596 

158.561 

-.176 

-.161 

106.685 

94 

.089 

486.433 

106.491 

158.565 

-.187 

-.175 

106.588 

95 

.101 

487.528 

106.396 

158.560 

-.201 

-.134 


9 6 

.114 

488 . 6 B 9 

106.308 

158.571 


-.202 

106,419 

97 

.128 

489.918 

106.232 

158.574 

-.237 

-.224 

106.349 

98 

.143 

491.214 

106.172 

158.576 

-. 248 

-.237 

106.295 

99 

.158 

492.574 

106.122 

156.578 

-.255 

-.243 

106.251 

ICO 

.174 

493.999 

106.06 6 

158.530 

-.269 

-.259 

106.220 

101 

.191 

495.487 

106.061 

158.532 

-.280 

-.269 

106.201 

102 

.209 

497.037 

106.049 

158.533 

-.287 

-.280 

1 06.193 

103 

.227 

498.649 


158.534 

-.2 96 

-.284 

106.202 


.246 

600.321 

106.060 

158.585 

-.314 


10 6.212 


.266 

502.052 

1 06.090 

158.535 

-.323 

-.321 

106.245 

106 

.236 

503 . 84 ? 

106.144 

158.535 

-.326 

-.321 

106.302 

1 07 

.307 

505.689 

106.194 

158.534 

— . 336 

-.329 

106.355 

108 

.329 

507.592 

106.271 

158.531 

jh 

|DQQ 

106.435 


.351 

509.549 

106.364 

158.578 



106.530 

110 

.374 

511.560 

106.476 

158.575 


I9Ebu 


111 

.397 

513.623 

106.611 

156.571 

-.368 

-.365 

106.731 

112 

.421 

515.736 

106.767 

158.566 

-.373 

-.374 

106,939 

113 

.446 

517.899 

106.955 

158.560 

373 

-.376 

107.129 

114 

.471 

520.111 

107.168 

158.553 

-.372 

-.376 

107.342 

115 

.496 

522.369 

107.421 

158.544 

-.359 

-.375 

107.595 

116 



iHrllrUl 


-.329 

-.350 

1 07 . P 94 


.549 

527.023 

108.060 

158.521 

-.288 

-.314 

10 H .231 

118 

.576 

529.415 

108.452 

158.506 

-.236 

-.269 

108.620 

119 
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Table 2 (Concluded) 


• 632 

534.322 

109.405 

156.465 

-. 095 

-.135 

109.566 

121 

• 660 

536.634 

109.941 

156.433 

-.020 

-.059 

lie . 098 

122 

• 689 

539.384 

110.515 

158.397 

.053 

.016 

110.668 

123 

• 719 

541.970 

111.099 

158.358 

.119 

.087 

111.246 

124 

• 748 

544.591 

111.686 

158.316 

.178 

.149 

111.830 

125 

• 778 

547.246 

112.259 

156.271 

.231 

• 205 

112.400 

126 

.909 

549.939 

11?. 813 

156.225 

• 280 

• 256 

112.951 

127 

.84 0 

552.663 

113.327 

158.180 

.323 

.304 

113.464 

126 

.871 

555.421 

113.777 

158.136 

.359 

.342 

113.914 

129 

.90 3 

558.215 

114,162 

158.110 

.387 

.376 

114.301 

150 

.935 

561.044 

114.433 

156.106 

.397 

.399 

114.577 

131 

.967 

563.906 

114.555 

158.119 

.361 

.396 

114 . 7 C 4 

152 

i.odo 

566.732 

114.461 

156.144 

.261 

.329 

114.672 

133 


Fig. 10. The airfoil sections have considerable openness. One engineering 
procedure to close the gap is to rotate the lower surface about the leading 
edge. The inverse procedure in its present form needs a more rigorous trailing 
edge closure model. Two such candidate procedures are described in the next 
section as recommendations for further study. 

RECOMMENDATIONS FOR FURTHER STUDY 
Trailing Edge Closure 

When a favorable pressure distribution is prescribed, it doesn's guarantee 
the resulting trailing edge thickness distribution to come out satisfactorily. 
To some extent, the trailing edge thickness can be controlled by adjusting the 
leading edge shape or the velocity potential value at the leading edge. 
Procedures to implement these ideas are described here. 

Leading Edge (Nose) Shape Alteration 

In a mixed analysis-inverse problem where the shape near the nose is 
usually prescribed and the objective is to weaken the shock or move the shock 
downstream, the shape of the nose can be used to control the trailing edge 
thickness. First, specify the nose shape given locally by y = Sqx”, where 
n and ag are two free parameters, and specify the desired Cp on the rest of 
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Fig. 6. Modified shape to weaken the shock system, n = 0.6932, M = 0.86, 
a = 4.86° 
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BASE WING GEOMETRY 

SPANWISE DIMENSION - INCHES 
40 80 120 160 200 240 280 320 360 



TEST STRIP 


Fig. 7. Base wing geometry with a test strip where 
wing modification is required 


BASE WING ROOT AIRFOIL 



Fig. 8. Base wing airfoil geometry 
at the root section 




DESIRED MIDSPAN PRESSURE DISTRIBUTION 



Fig. 9. Desired streamwise pressure distribution 
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X/C 


10. Modified airfoil shapes at two different span stations to 
provide laminar flow control pressure distribution at 
midspan 



the surface. The free parameters ag and n will then have to be adjusted using 
a gradient approach to satisfy a preset trailing edge thickness constraint. 

This method could possibly involve a mismatch in pressure at the point of 
transition from the analysis nose region to the Cp prescribed inverse region. 

If a mismatch in pressure occurs, then the specified pressure in the transition 
region must be allov/ed to vary to preserve a smooth pressure distribution. 

Leading Edge Velocity Potential Alteration 

(41 

In the case of a small disturbance methodology' the trailing edge 
closure was obtained by an alteration of the nose velocity potential. Such 
a procedure can also be tried in this full potential formulation. Let us define 
t|^ to be the trailing edge thickness for the k " span station. The objective 
then is to drive this tj^ to some preset value (will be zero for closed trailing 
edge airfoil) by perturbing the velocity potential ()> residing at the leading 
edge which is denoted by <(>mQ 5 £- Denoting the functional relationship of t|^ as 
^k “ ^k(*^N0SE)’ write the following expansion 


t 


k 




*i<('*'nose) * 


‘('*'nose)„ * ••• 

m 


m = kg, ••*, k-1, k, k+1, •••, k^ 
k = kg, •••, k-1, k, k+1, ••*, kg 





r (18) 


where kg and k^ denote the 
design mode, respectively, 
side of Eq. (18) is set to 


for all 


NOSE, 


m 


first inboard and final outboard span 
For trailing edge closure condition, 
zero which yields enough equations to 


stations under 
the left hand 
uniquely solve 




(19) 
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The Ml 


I'MnccI solution vector from Eq. (19) gives the amount of alteration 

c nuot ' ( j ( ) ( ) 

to be made on |4>(^0SE( drive |t|^| to zero. In Eq. (19), the | j symbol 
denotes a vector and [ ] ^ denotes the inverse of a matrix. Each element 
of this matrix is a partial derivative and a complete evaluation of all the 
matrix elements and the subsequent matrix inverse can be very costly and 
time-consuming, especially if several span stations are under design mode. 

To substantially reduce the computer time in the evaluation of matrix elements 
in Eq. (19), some tricks are used. First, the span station which has the 
maximum openness or fishtail is selected. For this span station (call it 


at,. 


k=k^) the influence function g 


(‘*’nose). 


is generated and that influence 


function distribution is kept the same for all other design span stations but 
the magnitude is scaled by the following 



( 20 ) 


It is recommended that both these procedures be tried in the currently 
developed inverse program. Besides the trailing edge closure model, further 
work is also recommended to assess the importance of pressure constraints in 
the inverse setting and the relationship between prescribed pressure and the 
freestream Mach number. 
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CONCLUSIONS 


An inverse procedure based on the full potential equation in conservation 
form has been developed for use in recontouring a given wing to produce a 
prescribed favorable pressure distribution. A density linearization scheme 
is introduced to aid in the application of the pressure boundary condition. 

The inverse logic is incorporated into the existing finite volume FL030 
analysis computer program. The new shape information is obtained from a mass 
flux integration procedure. The method is reasonably inexpensive and can be 
effectively used for shockless or shocked flow wing design. Tv/o procedures 
to control the trailing edge are proposed for further study. 
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APPENDIX A 


DERIVATION OF EQ. (11) IN THE DENSITY LINEARIZATION SECTION 

For simplicity, the derivation shown here is for two dimensions and the 
extension to three dimensions is straightforward. Considering only the 
directions, the contravariant velocities U and V can be written as 


U = + vCj, = A, 4.^ + 

V = Un^ + vriy = Aj*^ + A,*^ 

where 


(A-l) 


*1 = 'x ^ 

*2 ' 'x\ * '/y 

''3 = -lx ^ 

and u and v are the Cartesian velocity components along x and y. Using 
Eq. (A-l), the expression for density can be written as 

r vi i1/(y-i) 

P(4>) = [1 - ^ M2(A^<})2 + 2A24>^({>^ + A3<j)^ - 1)J . (A-2) 


The change in density due to small changes in the velocity potential cj) can 
be analyzed by substituting (<}) + 6(j)) for <{) in Eq. (A-2). 
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r Y-i 1 nV(Y-i) 

p(4>+ <S4>) = ~ 2 ~ + 2A2((j)^+ 6(f)^)(<{)^+ 6(j)^) + A2(<j)^+ 1 jj 

= ^ f^^|Ai<l>^ + 2A2(t)^(()^ + A3(()2 + 2A^(j)^6ct)^ + 2A2(<j)^54)^ + (|)^6({)^) 

|V(y-i) 


+ 2A3<f^6<j)^ - 1 IJ 
(y-l)M,2 


^ + 4>^<Sct)^) + A3c|)^6ct, J j 


-||1/(y-1) 


p((J)+6<J)) = p((|))|l - j 


(A-3) 


6p = p((}) + 6(l)) - p{(|)) = - P^"^M2 [u6({) +V6(J) I 

L ^ T| «J 


)2-Ym2 




64) 


differential operator 



(A-4) 


Equation (A-4) is the two dimensional analog of Eq. (11) in the main report. 
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APPENDIX B 


INSTRUCTIONS FOR THE USE OF THE INVERSE CODE IN ITS PRESENT FORM 

In subroutine XSWEEP (after XSWEEP.15), the user specifies the following 
information. 

a. IDU '>^1 index of the first upper surface point from the leading 

edge where inverse calculation starts. 

b. IDL I index of the first lov/er surface point where inverse 

calculation starts. 

c. KNIB First inboard span station index K where inverse calculation 

starts. 

d. KOUTB Last outboard span station where inverse calculation ends. 

Specification of modified pressures at half node points under the 
dimensional array name CPD(I,K), in the main program after MAIN. 80. 

The format of the output is shown in Tables 1 and 2 and explanations of 
columns CD to (J) are given on pages 12 and 13, 
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